Automated Synthesis of Dynamically Corrected Quantum Gates 
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Dynamically corrected gates are extended to non-Markovian open quantum systems where lim- 
itations on the available controls and/or the presence of control noise make existing analytical 
approaches unfeasible. A computational framework for the synthesis of dynamically corrected gates 
is formalized that allows sensitivity against non-Markovian decoherence and control errors to be 
perturbatively minimized via numerical search, resulting in robust gate implementations. Explicit 
sequences for achieving universal high-fidelity control in a singlet-triplet spin qubit subject to real- 
istic system and control constraint are provided, which simultaneously cancel to the leading order 
the dephasing due to non-Markovian nuclear-bath dynamics and voltage noise affecting the control 
fields. Substantially improved gate fidelities are predicted for current laboratory devices. 

PACS numbers: 03.67.Pp, 03.67.Lx, 73.21. La, 85.75.-d 



I. INTRODUCTION 

Achieving high-precision control over quantum dynam- 
ics in the presence of decoherence and operational errors 
is a fundamental goal across coherence-enabled quantum 
sciences and technologies. In particular, realizing a uni- 
versal set of quantum gates with sufficiently low error 
rate is a prerequisite for fault-tolerant quantum compu- 
tation [1]. Open- loop control based on time-dependent 
modulation of the system dynamics has been extensively 
explored as a physical-layer error-control strategy to meet 
this challenge. Two main approaches have been pur- 
sued to date: on the one hand, if the underlying open- 
system relaxation dynamics is fully known, powerful vari- 
ational techniques and/or numerical algorithms from op- 
timal quantum control theory (OCT) may be invoked 
to optimize the target gate fidelity, see e.g. [2-5] for 
representative contributions. On the other hand, dy- 
namically corrected gates (DCGs) [(>] have been intro- 
duced having maximum design simplicity and portabil- 
ity in mind: close in spirit to well-established dynam- 
ical decoupling techniques for quantum state preserva- 
tion in non-Markovian environments [7], DCG sequences 
can achieve a substantially smaller net decoherence er- 
ror than individual "primitive" gates by making minimal 
reference to the details of the system and control specifi- 
cations. In principle, the use of recursive control design 
makes it possible for the final accuracy to be solely lim- 
ited by the shortest achievable control time scale [8] . Re- 
markably, simple DCG constructions underly the fidelity 
improvement reported for spin-motional entangling gates 
in recent trapped- ion experiments [9]. 

While obtaining a detailed quantitative characteriza- 
tion of the noise mechanisms to overcome is imperative 
to guarantee truly optimal control performance, this re- 
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mains practically challenging for many open quantum 
systems of interest. In addition, current approaches for 
applying OCT methods to non-Markovian environments 
rely on obtaining suitable simplifications of the open- 
system equations of motion (e.g., via identification of a 
finite-dimensional Markovian embedding [2] or approxi- 
mation through time- local coupled linear equations [5]) - 
which may be technically challenging and/or involve non- 
generic assumptions. Since in DCG schemes the error 
cancellation is engineered at the level of the full system- 
plus-environment Hamiltonian evolution, two significant 
advantages arise: environment operators may be treated 
symbolically, avoiding the need for an explicit equation 
of motion for the reduced dynamics to be derived; in 
contrast to error-control approaches designed in terms of 
gate propagators (notably, fully compensating compos- 
ite pulses for systematic control errors [10, 11]), work- 
ing at the Hamiltonian level allows to more directly re- 
late to physical error mechanisms and operational con- 
straints. Despite incorporating realistic requirements 
of finite maximum control rates and amplitudes, ana- 
lytic DCG constructions nonetheless rely on the assump- 
tion that complete control over the target system can 
be afforded through a universal set of stretchable con- 
trol Hamiltonians [(>]. This requirement is typically too 
strong for laboratory settings where only a limited set of 
control Hamiltonians can be turned on/off with sufficient 
precision and speed, and universality also relies on in- 
ternal always-on Hamiltonians. Furthermore, portability 
comes at the expenses of longer sequence durations, mak- 
ing DCGs more vulnerable to uncompensated Markovian 
decoherence mechanisms. 

In this work, we introduce a control methodology that 
results in an automated recipe for synthesizing DCGs via 
numerical search. This is accomplished by relaxing the 
portability requirement and utilizing the full details of 
the control. While the resulting "automated DCGs" (aD- 
CGs) are still synthesized without quantitative knowledge 
of the underlying error sources, they overcome the re- 
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strictive assumptions of analytical schemes and lead to 
drastically shorter sequences. As an additional key ad- 
vantage, our Hamiltonian-engineering formulation lends 
itself naturally to incorporating robustness against mul- 
tiple error sources, that can enter the controlled open- 
system Hamiltonian in either additive or multiplicative 
fashion. This allows for aDCGs to simultaneously cancel 
non-Markovian decoherence and control errors, as long 
as the combined effects remain perturbatively small. 

We quantitatively demonstrate these advantages by fo- 
cusing on a highly constrained control scenario - the two- 
electron singlet-triplet (S-To) spin qubit in GaAs quan- 
tum dots (QDs)[f2]. In spite of ground-breaking exper- 
imental advances [13-16], boosting single-qubit gate fi- 
delities is imperative for further progress towards scal- 
able quantum computation and is attracting intense the- 
oretical effort [17-19]. Recently introduced SUPCODE 
composite-pulse sequences [19], for instance, are (analyt- 
ically) designed to achieve insensitivity against decoher- 
ence induced by coupling to the surrounding nuclear-spin 
bath, however they do not incorporate robustness against 
voltage noise, which is an important limitation in exper- 
iments [16]. Here, we provide explicit aDCG sequences 
for high-fidelity universal control in S-To qubits, which 
cancel the dominant decoherence and exchange-control 
errors, while respecting the stringent timing and pulsing 
constraints of realistic S-To devices. The resulting se- 
quences use a very small number of control variables and 
a fixed base pulse profile, which streamlines their exper- 
imental implementation. Up to two orders of magnitude 
improvement in gate fidelities are predicted for parameter 
regimes appropriate for current experimental conditions. 



II. CONTROL-THEORETIC SETTING 

We consider in general a d-dimensional target quantum 
system S coupled to an environment (bath) B, whose 
total Hamiltonian H on Hs ® Hb reads 



a way that, without loss of generality, we may express 



H = [H s + H ctlX0 (t)]®I B +H e 

H e = H e ,int + He : ctrl(t), 



(1) 



where I B (s) denotes the identity operator on B(S), 
Hs accounts for the internal ("drift") system's evolu- 
tion in the absence of control, and the time-dependent 
H c tii,o(t) represents the intended control Hamiltonian on 
S. The total "error Hamiltonian" H e encompasses the 
bath Hamiltonian, unwanted interactions with the bath, 
as well as deviations of the applied control Hamiltonian 
from Hctr^oWi subject to the requirement that the un- 
derlying correlation times are sufficiently long. Formally, 
we require that \\H e \\ <§; ||ff(i)||, where ||X|| is the oper- 
ator norm of X = X^ (maximum absolute eigenvalue of 
X) [6, 7]. In order to "mark" the error sources, we char- 
acterize the strength of each independent contribution to 
-f/^int in terms of dimensionless parameters {S a }, in such 
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Hb + H SB — 22 3aS a 



with S a being a Hermitian operator basis on Hs and 
B a acting on Hb, respectively, and the bath internal 
Hamiltonian Hb = SqBo- We assume that the B a are 
norm-bounded but otherwise quantitatively unspecified. 
In particular, if B a are treated as scalars (B a = ioJ-B), 
we may formally recover the limit of a classical bath, 
whereby Bq — and the system Hamiltonian is effec- 
tively modified in a random (yet slowly time-dependent) 
fashion. Note that, as long as we are interested in cancel- 
ing effects that are first order in the error sources, there is 
no distinction between the B a being actual operators or 
scalars. Similarly, we characterize the independent error 
sources in H ctr ifl(t) by letting 

#e,ctrl(i) = ^SpHp (t)<8 Is, 

where Hp(t) are known system operators, while the pa- 
rameters 6p remain unspecified. For notational conve- 
nience, we shall label all the unknown parameters sym- 
bolically and collectively by 5 = {S a , 8p}. 

In an ideal error-free scenario, 5 = 0, the system 
evolves directly under the action of the control, in the 
presence of its internal drift Hamiltonian. We assume 
that in this limit, S is completely controllable, that is, ar- 
bitrary unitary transformations on S can be synthesized 
as "primitive gates" by suitably designing ifctri,o(0 m 
conjunction with Hs- As mentioned, we are particularly 
interested in the situation where the latter is essential for 
controllability to be achieved [20]. The available control 
resources may be specified by describing 

flrtrl,o(*) =J2 Ca ^ Ha ® lB > 

a 

in terms of the admissible (nominal) control inputs and 
Hamiltonians. Beside restrictions on the set of tun- 
able Hamiltonians H a , limited "pulse-shaping" capabil- 
ities will typically constrain the control inputs c a {t) as 
system-dependent features of the control hardware. For 
concreteness, we assume here that £/ctri,o(i) is decom- 
posed as a sequence of shape-constrained pulses applied 
back to back and also constrain pulse amplitudes {hi} 
and durations {r^} to technological limitations such as 

^min — hi 5~ ^max, 7"miii — T% — ^"max- 

Ideally, if the target unitary gate is Q, the objective for 
gate synthesis is to devise a control Hamiltonian H ctI \ } o(t) 
such that (up to a phase), 



Q = Tcxp 



i I (H s + Hcttifl(s))ds 



(2) 



where T denotes time ordering and r is the running time 
of the control. The ideal evolution naturally defines a 
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toggling-frame unitary propagator given by 



UQ(t) = Texp 



i I (H s + H ctT i t0 (s)) ds 
o 



(3) 



which traces a path from Is to Q over r. If H e ^ 0, ap- 
plication of i? c tri,o(i) over the same time interval results 
in a total propagator of the form 



U. 



(*) = jj(0) 



Q[r] 



■&]«P(-«E$ T] ). 



where Uq)-, = Q and Eq), is an "error action" operator 
on Us ®%b that isolates the effects of undesired terms 
in the evolution [(>]: 

exp(-iE^ r] ) = Texp U Q {sfH e U Q {s) ds] . (4) 

The norm of the error action can be taken to quantify 
the error amplitude per gate (EPG) in the presence of 5. 
The EPG in turn upper-bounds the fidelity loss between 
the ideal and actual evolution on S once its "pure-bath" 
components are removed. More concretely, define 

mod B E = E-I S ® Tr s (E)/d, 

that is, a projector that removes the pure-bath terms in 
E (note that vho&bE = E if E is a pure-system operator 
of the form A(g)Ig, as for a classical bath). Then the fol- 
lowing (not tight) upper bound for the (Uhlman) fidelity 
loss holds independently of the initial states [8, 21]: 



l-fu< \\mod B E i 



(<5) 



Thus, reducing the EPG can be used as a proxy for re- 
ducing gate fidelity loss. While for a primitive gate im- 
plementation the EPG scales linearly with S, the goal of 
DCG synthesis is to perturbatively cancel the dependence 
on S in ^gr T i up to a desired order of accuracy, to realize 
the gate in a manner that is as error free as possible as 
long as 6 is small. For simplicity and immediate appli- 
cation, we focus here on first-order aDCG constructions, 
for which \\mod B E Q s) [r] \\ = 0(S 2 ). 



III. SYNTHESIZING DYNAMICALLY 
CORRECTED GATES 

A. Existential approach 

Recall that two main requirements are required in first- 
order analytical DCG constructions [6, 8]: (i) primitive 
gate implementations of the generators of a "decoupling 
group" associated with the algebraic structure of EPGs 
and (ii) particular implementations of the target gate 
Q (as Q*) and the identity gate (as Iq) as sequences 
of primitive gates such that Q* and Iq share the same 
first-order EPG, making them a "balance pair". While 



(i) is provided by controllability and leads directly to a 
constructive procedure for correcting to the first-order 
the identity evolution, (ii) is essential for modifying this 
procedure in such a way that the net first-order error 
cancellation is maintained, but Q is effected instead. 

Generating balance pairs require further adjustment 
of gate control parameters to form a controllable rela- 
tionship between EPGs of gate implementations, holding 
as an identity regardless of the value of 6 (or B a ). For 
example, in the absence of drift dynamics and control er- 
rors, such a controllable relationship can be engineered 
by "stretching" pulse profiles in time while the amplitudes 
are reduced proportionally, resulting in different realiza- 
tions of the same target, with EPGs that scale linearly 
with the gate duration. Similarly, in the presence of a 
multiplicative control error, primitive gates with phys- 
ically equivalent (modulo 2tt) angles of rotation result 
in different EPGs (note that similar geometric ideas are 
used in designing composite pulses [11]). We argue next 
that knowing the control description and marking the er- 
ror sources {5j} does still lead to (ii) as long as control 
constraints allow us to tap into a continuum of different 
gate implementations. 

The multitude of pathways for realizing a primitive 
gate Q increases with gate duration/subsegments as a re- 
sult of availability of more control choices and ultimately 
a simpler control landscape [22]. Assume that (Al) such 
primitive implementations may be parametrized as Q[t\. 
We aim to show that a balance pair or, alternatively, 
a direct cancellation of the EPG of Q may be found. 
The gist of our argument is most easily given for a single 
qubit, with the Pauli operators chosen as the operator 
basis {S a } for error expansion. Using the fact that the 
interactions among different error sources can be ignored 
up to the first order, the basic idea is to start with the 
first error source, Si, and then use the resulting gates re- 
cursively for the next error source until all error sources 
are exhausted. Ignoring error sources other than Si, let 
us thus expand E Q 5 ^ = Si J2 a eQ[r], a S a ® B x . 

Assume in addition that (A2), as a function of the pa- 
rameter r, the range of the real- valued functions e.Q\ T ], a 
extends to infinity in positive or negative directions. Con- 
sider now "projection blocks" composed of two Pauli gates 
applied back to back, that is, 

-Pab~a] = $a \^~a ] $a \j~ot ] , 

with a corresponding EPG given by 2(5ies Q [ Tci ] jCt >S , a ® Bi, 
which is purely along S a . By virtue of (a2), we can find a 
continuum of (r, T a ) pairs such that 2eg r T i Q = ±eQr T i a 
for all Pauli directions a, meaning that we may reproduce 
each error component in Eq^ up to a sign. Those Pauli 
components a_ that reproduce error with a negative sign 
are combined in sequence with Q[t] to form a longer gate 



Q* = Q[t] l[P a [r a 



If all Pauli components can be matched with negative 
signs, the resulting gate will cancel all error components 
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and a DCG construction is provided by Q* . Otherwise, 
the Pauli components a + that are matched with positive 
sign, are combined to produce an identity gate, 

a + 

that matches the error of Q* . Hence, (Q*,Iq) form a 
balance pair and can be used to produce a continuum of 
constructions of a DCG gate Q^> [r] that cancel the error 
source 8%. Provided that the assumption (A2) remains 
valid for this new composite constructions, we can repeat 
the procedure to remove the other error sources. 

We remark that assumption (A2) essentially implies 
that the domain of the errors as a function of implemen- 
tation parameters for a fixed unitary gate is not compact, 
so that arbitrary magnitudes of each error component can 
be sampled by appropriately choosing the implementa- 
tion parameters. Such arbitrary large domains need not 
not exist in the primitive gate implementations (natu- 
rally or due to constraints), or only discrete error values 
may be reachable. Nonetheless, we may still enlarge the 
accessible range of errors for the target gate Q by attach- 
ing a continuously parametrized family of identity gates. 
Universal controllability of the system implies that not 
only any gate U but also its inverse U~ x may be reached. 
Implementing U, followed by its inverse U -1 , produces an 
implementation of the identity Ijj "parametrized" by the 
original gate U. In the absence of degeneracies (relation- 
ships between the errors that could be used separately 
to provide a balance pair), the EPG associated with Ejj 
has then a continuos domain. Clearly, applying Ijj fol- 
lowed by the target gate Q still realizes the gate Q but 
the resulting EPG is now given by Eq + Ej v , which is 
parametrized by U. By applying sufficiently many copies 
of Ijj before applying Q, the error can be extended to ar- 
bitrary large domains as desired. 

B. aDCGs: Computational Approach 

The aDCG sequences generated by following the above 
existential argument tend to be far too long and com- 
plex to be useful in realistic control scenarios. Also note 
that in principle, the construction of single- or two-qubit 
DCGs in n-qubit registers may be handled similarly by 
using multi-qubit Pauli operators as a basis for EPG ex- 
pansion and for building projection blocks. However, the 
sequence complexity tends in this case to also grow ex- 
ponentially with n [8], making the need for more effi- 
cient synthesis procedures even more essential. Just as 
complete controllability provides an existential founda- 
tion to numerical OCT approaches for unitary gate syn- 
thesis when H e = [18], our argument legitimates a nu- 
merical search for aDCGs in the presence of H e . Since 
the objectives of gate realization and perturbative error 
cancellation are not inherently competing, the numerical 
search can be described as multi- objective minimization 
problem, as we detail next. 



Let the nominal control Hamiltonian i? c tri,o(0 be 
parametrized in terms of control variables {xi} and define 
objective functions as follows: 

F({x<}) = dist(C/g= 0) ,Q), (5) 

BjGjiixi}) = ||9mod B ^ T] /^|| tf=0) (6) 

where j labels independent error sources and Bj symbol- 
ically denotes bath operators that mark error sources in 
-ffe,int (recall that Bj = 1b for control error sources) to 
ensure that F and Gj only depend on the known quan- 
tities Xi. Minimizing only the first objective, F = 0, 
corresponds to achieving exact ideal primitive gate syn- 
thesis, Eq. (2). As an appropriate distance measure for 
unitary operators in Eq. (5), we use 

disb(U,V) = [1 - \Tr{U*V)\/d\ 1/2 , (7) 

which is a standard phase-invariant choice [18]. Minimiz- 
ing the objectives in Eq. (6) corresponds to first-order 
sensitivity minimization. Thus, solving for F = = 
Gj,Vj, results in an implementation of Q that is insen- 
sitive to the perturbative parameters Sj, yielding a robust 
control solution as long as Sj is small. 

Evaluating Gj apparently requires solving the full 
time-dependent system-plus-bath Schrodinger equation 
parametrized by the controllable pulse shapes. In fact, 
once the error sources (including bath operators) are 
treated as first-order symbolic variables, Gj can be eval- 
uated by effectively solving the Schrodinger equation on 
the system only, in order to determine the appropriate 
toggling-frame propagator, Eq. (3), and then evaluate 
the required error action, Eq. (4), by invoking a Mag- 
nus expansion [6]. Specifically, if the control variables 
Xi = {("Tj, hi)}, the sequence propagator reads 

where r = X)i T * anc ^ Ux^\s) is the i-th pulse propaga- 
tor corresponding to the variable x% at time s, with its 
associated first-order error action, 

4?.= rug)(s)iH e Ug>(s)ds. 
Jo 

To the first order in 6, the total EPG is in turn given by 

i 

where Vi = ■ • ■ t/i? denote the "partial" product 

of gate propagators up to and excluding the i-the gate 
[6]. While, as noted, for a first-order aDCG the resulting 
accuracy modgE^^ \ \ = 0(S 2 ), a more quantitative esti- 
mate of the actual conditions of applicability requires es- 
timating the dominant uncorrected second-order errors. 
Technically, this can be carried out by means of standard 
algebraic techniques however is not straightforward [2'->] 
and beyond our present scope. Instead, we focus in what 
follows on addressing the construction and performance 
of first-order aDCGs in concrete illustrative settings. 
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IV. APPLICATION TO SINGLET-TRIPLET 
QUBITS 

Consider first the following single-qubit specialization 
ofEq. (f): 



Ba x + J (t)a z 



(8) 



H e = I s ®5 B + a x ®8 1 B x + J (t)8 2 a z 



where the operator- valued 5±B X and the system drift B 
couple to the system along a x and the nominal control 
Jo(t) and a multiplicative error 82 couple along a z . Al- 
though explicitly included, the bath internal Hamiltonian 
does not play a role in the first-order removal of decoher- 
ence and is automatically accounted for in mods. On 
the other hand, the drift term B is essential for complete 
controllability and analytical DCG constructions are not 
viable even in the limit 82 — > 0. Thus, the need to ef- 
fectively address both noise sources Si , 62 for a generic 
operating point B mandates the use of aDCGs. 

While useful as a template for single-axis control sce- 
narios in the presence of internal drift and dephasing, a 
semi-classical version the above model is relevant, in par- 
ticular, to describe a universally controllable S-To qubit. 
In this case, the logical qubit subspace is spanned by 
{\S), To)}, the singlet and triplet states of two electrons 
on separate QDs [14, 15, 19, 24] and, provided that the 
number of bath nuclear spins is sufficiently large [24] , the 
following simpler Hamiltonian is appropriate and widely 
used for this system [18, 19]: 



H 



B + SB(t) 



1 



J{t)a z 



(9) 



Physically, the drift term B is a known static magnetic 
field gradient between the two QDs that includes an 
Overhauser field from the nuclear spin bath, SB(t) (cor- 
responding to 2SiB x ) accounts for random fluctuations 
of B due to coupling to nuclear flip-flop processes [25], 
and J(t) is the exchange splitting. In practice, J(t) is 
tuned by control of an electrostatic gate voltage [14], and 
voltage fluctuations due to charge noise result in a noisy 
control Hamiltonian of the form J(t) — Jo(f)(l + SJ(t)), 
where 5J(t) thus corresponds to 282 ■ We assume that 
both noise sources may be treated as Gaussian quasi- 
static processes, with their "run-to-run" distribution be- 
ing characterized by standard deviations a$B and <j$j. 
While in practice the noise is not completely static, we 
expect our considerations to remain valid as long as high- 
frequency noise components decay sufficiently fast and 
the resulting aDCGs are short relative to time scales 
over which white charge noise may become important. 
Phenomenologically, the dephasing induced by the fluc- 
tuating Overhauser field is consistent with a power-law 
noise spectrum of the form S(w) ~ lo~ 2 over a wide 
spectral range [26]. Likewise, recent experiments indi- 
cate that voltage noise also arises overwhelmingly due to 
low- frequency components with an approximate 1/ f de- 



cay at low operating temperatures [16]. From experimen- 
tally measured values of TJ , we use here o~sb/{2tt) ^ 0.15 
MHz [13] and a SJ < 1/50 [27]. 

In constructing aDCGs, we shall choose values of the 
internal drift (B) and of the nominal control field (Jo(£)) 
appropriate for the QD setting of Eq. (9). We stress, 
however, that the same solution is found from (and ap- 
plies to) the fully quantum model Hamiltonian of Eq. 
(8). In practice, the drift term B can be set to a fixed 
value B/2n e [0.03, 0.2] GHz, which we choose at 0.1 GHz 
[28]. The control field Jo(t)/2ir is taken to be positive 
and smaller than J max = 0.3 GHz. We recognize the 
finite rise, delay, and drop times associated with pulse 
generators by fixing a pulse profile. Thus, during each 
pulse, with time t' measured from the pulse start, the 
control function Jo{t') is given by hic(t'/ri), where c(x) 
is the pulse shape function. We digitize the pulse shape 
function for numerical evaluation. In contrast to merely 
bounding the pulse times and allowing pulse durations 
as extra control variables, we enforce the pulse durations 
to be fixed at r\ = r = 3 ns, that is compatible with 
the currently most widespread pulse generators temporal 
resolution of 0.83 ns. The search space of the pulse am- 
plitude control variables is thus given by Xi = hi. While 
removing Tj from the control variables results in more se- 
vere constraints, it also corresponds to a reduction of the 
search space. We verified that all of our results were re- 
produced with variable but lower-bounded pulse widths 
as well. The objective functions F, G\ and G2 are com- 
puted explicitly in terms of each constituting pulse pa- 
rameter hi according to the general procedure described 
in Sec. III.B. 

For the resulting multi-objective minimization, we in- 
troduce numerical weight factors Ai and A 2 and form a 
single objective function 

0({hi}) = F + Aid + A 2 G 2 . 

Choosing small values of Ai (= 10~ 5 ) work best in direct- 
ing the search from solutions that synthesized the target 
gate only at (F = 0) first, and then towards the error- 
corrected solution F = 0, G\ = = G 2 . Motivated by 
our existential argument, the intuition is to avoid the 
local minima associated with multiple objectives and fo- 
cus on a single objective which, once realized, will give 
weight to the other objectives iteratively. We solve each 
aDCG search problem using off-the-shelf (Matlab's fmin- 
CON function) search routine (within minutes), with the 
default choice for solving constrained nonlinear optimiza- 
tion without specifying a precalculated gradient or Hes- 
sian. We start the search with a small number of pulses, 
n, which is then incremented until the minimal value of 
the objective function comes close to the machine pre- 
cision (« 10~ 16 ). Fig. 1 (top) depicts the synthesized 
control profiles for a universal set of single-qubit aDCGs. 

Once aDCG sequences are found, evaluating their ef- 
fectiveness for the S-T qubit can take direct advantage of 
the effectively closed-system nature of the model Hamil- 
tonian in Eq. (9) , thus avoiding the need of explicit spin- 
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FIG. 1: (Color online) Upper panel: aDCG control profiles 
for exp(— ina z /8) (top) and the Hadamard gate (bottom). 
Lower panel: Fidelity loss for uncorrected vs. corrected 
gates, evaluated as I — / = dist((7g 5S '' 5J ', Q) 2 [Eq. (7)], for 
Q = exp(— inaz/8) as a function of relative magnetic field 
gradient error SB/B (S/(2tt) = 0.1 GHz). SJ = corre- 
sponds to perfect exchange control, whilst SJ = 0.01 is close 
to current realistic levels. The fact that SJ = 0.01 is fixed 
is responsible for the eventual performance plateaux where 
the latter dominates over SB effects. Nevertheless, the aDCG 
advantage is always maintained. 
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(n = 17 pulses). The higher slope of the fidelity loss as 
a function of SB when S J = is the signature of a per- 
turbative error cancellation and the aDCG advantage is 
maintained even with SJ > 0, implying robustness with 
respect to both error types. 

In order to make contact with experimentally relevant 
ensemble-averaged fidelities, we further evaluate the av- 
erage of single-run fidelities f(SB, SJ) with respect to 
noise realizations, by assuming that SB and SJ are inde- 
pendent and normally distributed random variables with 
variance o~\ B and o~\j- The results are summarized in 
Fig. 2. Both noise sources adversely impact the expected 
gate fidelity but aDCGs are far less affected, resulting in 
robust gates roughly as long as ctsb/b + o~8J ~ 0.1. 



V. CONCLUSION 

Our procedure can be interpreted as a automated gate 
compiler which incorporates detailed information about 
the controllable parameters and their range of opera- 
tions, along with qualitative information about the error 
sources affecting the evolution. Compared to mere (prim- 
itive) gate synthesis, the resulting increase in complexity 
scales proportionally to the number of error sources. Our 
approach applies to any Hamiltonian control setting, and 
for weak enough error sources, even higher-order cancel- 
lation can be achieved in principle. 

Thanks to the slow dynamics of the nuclear spin bath 
and fast control pulses available, electron spin qubits pro- 
vide an ideal experimental testbed for validating our ap- 
proach. While additional experimental details may be 
captured in more sophisticated ways, we believe that our 
framework is general and flexible enough for its effec- 
tiveness not to be compromised. In particular, further 
analysis is needed to quantify the effect of white elec- 
trical noise on aDCG sequences, as well as to possibly 
minimize its influence by penalizing large values of the 
exchange splitting in the numerical search. It is thus our 
hope that significantly improved single-gate fidelities will 
be achievable in S-To qubits by aDCG sequences that op- 
erate under realistic noise levels and control limitations. 



FIG. 2: (Color online) Ensemble averaged fidelity loss as 
a function of relative magnetic field-gradient and exchange- 
control noise for uncorrected (top, dark) vs. corrected (bot- 
tom, light) implementation of Q = exp(— i-Ka z /8). The solid 
(black) lines on each surface correspond to typical values 
osj ~ 0.02 and ctsb/b ~ 10 -3 for current S-To devices. 

bath simulations and quantum process tomography. Fig. 
1 (bottom) depicts the fidelity loss for an uncorrected 
(n = 3 pulses, obtained through the same numerical pro- 
cedure with Ai = A 2 = 0) vs. corrected implementation 
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